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The finite temperature transition of QCD with two degenerate light quarks 
is studied on the lattice with a renormalization group improved gauge action 
and the Wilson quark action. We have made simulations on an 8 3 x 4 lattice 
near the chiral transition point. It is shown that the chiral condensate which 
is the order parameter of chiral symmetry satisfies remarkably a scaling re- 
lation with the exponents of the three dimensional 0(4) Heisenberg model. 
This indicates that the chiral transition in two-flavor QCD is of second order 
in the continuum limit. 



Understanding the nature of the finite temperature transition in QCD 
with two degenerate light quarks (Np = 2) is an important step toward 
the clarification of the transition in the real world. It is plausible from an 
argument based on an effective o model [|IJ that if the chiral transition (tran- 
sition in the chiral limit) is of second order, QCD with two flavors belongs to 
the same universality class as the three dimensional 0(4) Heisenberg model. 
This universality provides us with several useful scaling relations that can be 
confronted with numerical results. These scaling relations were first tested 
on the lattice for staggered quarks || , where some evidences consistent with 
the 0(4) scaling were reported. 

In this paper, we study scaling behavior of the QCD transition with two 
degenerate light quarks on the lattice using the Wilson formalism of fermions, 
and compare the result with a scaling relation conjectured. The Wilson for- 
malism of fermions is the only known lattice formalism which possesses a 
local action for any number of flavors. Therefore, it is important to investi- 
gate the scaling behavior with Wilson quarks and compare the results with 
those for staggered quarks for which the action is non-local for two flavors. 

We use a renormalization group improved gauge action || which is ex- 
pected to be closer to the renormalized trajectory compared with the stan- 
dard one-plaquette gauge action, and therefore is expected to show scaling 
behavior on a coarse lattice. 

The continuum limit does not depend on the choice of the action if the ac- 
tion belongs to the same universality class. However, the way of the approach 
to the continuum limit does depend on the choice of the action. When one 
uses the standard one-plaquette gauge action and the Wilson quark action 
for a study of QCD with two flavors at finite temperatures, one encounters 
the existence of severe lattice artifacts on lattices with the lattice size in the 
euclidean time direction N t = 4 and 6 [|, ||, |6| . Note the temperature on the 
lattice is given by T = l/(aN t ) with a being the lattice spacing. A way out 
of these lattice artifacts with the standard action is to increase N t in order to 
further approach toward the continuum limit. However, our previous study 
on an 18 3 x 24 lattice suggests that Nt > 18 is required to have the chiral 
transition in the scaling region 0. It is difficult to perform simulations on 
such a large lattice even with a powerful computer of today. Therefore, we 
instead adopt a renormalization group improved gauge action and the Wilson 
quark action. 

Before going into the discussion of the present work, let us first describe 
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unexpected phenomena observed by the MILC collaboration with the stan- 
dard gauge action. The transition is smooth both for heavy and light quarks, 
while the transition is very sharp at N t = 4 (even a first order transition at 
N t = 6) in the range of intermediate quark masses (J3 ~ 5.0) jp. This is 
completely different from what is supposed to be realized in the continuum 
limit: as the quark mass increases from the chiral limit, the transition be- 
comes weaker and it becomes strong again only when the quark mass is heavy 
enough to recover the first order transition of the SU(3) gauge theory. They 
also observed a cusp in the behavior of m 2 in terms of 1/K with {3 fixed (m n 
is the screening pion mass and K is the hopping parameter). Furthermore, 
at f3 < 5.3, the value of the quark mass m q , defined through an axial- vector 
Ward identity || |9|, depends on the phase of the system and shows unex- 
pected 1/K dependence with (3 fixed in the high temperature phase |4], [II] 



which is in sharp contrast with the case of (3 > 5.5 where it does not de- 
pend on the phase (Tl|. In the following, we show that these unexpected 



phenomena are removed by improving the lattice action. 

Now let us discuss the present work. The renormalization group improved 
gauge action S™ we use is given by 

S T = ~2 ( c o EC 1 >< 1 lo °P) + c i EC 1 >< 2 lo °P)} (!) 

with c\ = —0.331 and Co = 1 — 8ci ||. Here the loops in the sums are 
defined by the trace of ordered product of link variables, and each oriented 
loop appears once in the sum. 

We mainly perform simulations for N t = 4 with spatial lattice size 8 3 at 
(3 = Q/g 2 = 1.0 — 4.0 near the finite temperature transition point. We use 
an anti-periodic boundary condition in the t direction and periodic boundary 
conditions otherwise. The configurations are updated using the hybrid Monte 
Carlo method with the molecular dynamics time step size 5t = 0.01 (except 
for a simulation point (3 = 1.5, K = 0.206 on the N t = 4 lattice where 
5t = 0.002). We perform simulations of about 100 — 1300 trajectories after 
thermalization for each set of parameters, (3 and K. Errors are estimated by 



the jackknife method. A preliminary report is given in Ref.[|l(| . 

Our results for the phase diagram are shown in Fig. |l]. The chiral line K c 
U is defined by the vanishing point of the screening pion mass, = 0, on 
an 8 4 lattice with periodic boundary conditions (doubled in the t direction 
for the spectrum calculation). The finite temperature transition/crossover 
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line K t for N t = 4 is determined on the 8 3 x 4 lattice from the condition 
that P = 0.10(2) where P is the Polyakov line expectation value. This 
condition corresponds to the criterion that the crossover point is identified 
as the peak position of the susceptibility of the Polyakov loop, although our 
limited statistics sometimes makes the peak of the susceptibility not so clear. 
We also note that this is consistent with the criterion that the crossover 
points are identified from rapid changes of physical quantities (see Fig. |2|). 
We identify the crossing point of the K c and K t lines as the point of the finite 
temperature chiral transition point, (3 ct (cf. Ref. 0). 

For the case of the standard gauge action, we pointed out that the simu- 
lation points where the strong transitions mentioned earlier are observed on 
the K t line || are just those where the K t line approaches toward the K c line 
after initially deviating from it due to the cross-over phenomenon between 
weak and strong coupling regions of QCD, and that therefore it seems plau- 
sible that the strong transition is a result of lattice artifacts caused by this 
unusual relation of the K t and K c lines f| ||, [H| . 

Unlike for the case of the standard gauge action, the distance between 
the K t and K c lines shown in Fig. [I] grows monotonically when we increase 
(3 from the chiral transition point (3 c t ~ 1.4. In accord with this, the tran- 
sition becomes monotonically weaker with (3 as shown in Fig. ^|a for the 
Polyakov loop. This is in sharp contrast with the case of the standard gauge 
action. Note that the change of the Polyakov loop at the finite temperature 
transition/crossover is very smooth for a wide range of (3. The change of 
is also smooth for larger values of (3 and it becomes sharper as (3 de- 
creases (cf. Fig. |||b). The straight line envelop of in the high temperature 
phase (Nt = 4) agrees with in the low temperature phase (Nt = 8), and 
corresponds to the PCAC relation oc m q . The smoothness of physical 
observables strongly suggests that the transition is a crossover at (3 > (3 ct . 

It should be noted that the value of rri^ on the K c line monotonically 
decreases to zero as (3 — ► (3 ct + (cf. Fig. |^), suggesting that the chiral 
transition is continuous. We also note that on the K t line shows a similar 
monotonic decrease [TIJ. It might be emphasized that the value of m q in the 
high temperature phase agree well with that in the low temperature phase 
and does not show the strange behavior mentioned before for the case of the 
standard action (see Fig. |]). 

These nice properties which are in accordance with naive expectations 
encourage us to begin a scaling study with Wilson quarks. 
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From the universality argument we expect that magnetization M near 
the second order phase transition point can be described by a single scaling 
function: 

M/h 1 ' 5 = f(t/hVM) (2) 

where h is the external magnetic field and t = [T — T c (h = 0)]/T c (h = 0) is 
the reduced temperature. For three dimensional 0(4) models, the critical 
exponents in (|) are given by 1/(J3S) = 0.537(7) and 1/5 = 0.2061(9) 0. 
In QCD, h corresponds to the quark mass and M corresponds to the chiral 
condensate. DeTar tested this scaling for the case of two flavor staggered 
quarks [jn| and reported that data are consistent with 0(4) and 0(2) scaling. 

For Wilson quarks, the naive definition of for the chiral condensate 
is not adequate because the chiral symmetry is explicitly broken due to the 
presence of the Wilson term. A proper subtraction and a renormalization 
are required to obtain the correct continuum limit. A properly subtracted 
(tyty) can be defined via an axial Ward identity ||: 

(M) sub = 2m q aZj2(^M0)} (3) 

X 

where Z is the renormalization coefficient. This definition is consistent with 
the identification of the magnetization in Ref.]l|]. ( 1 I , \E') su b was shown to have 
a non-vanishing value in the chiral limit in the confining phase of quenched 
QCD For our purpose, it is enough to use the tree value: Z = (2K) 2 . 



Our results of ( 1 I r \[ / ) su b for N t = 4 are shown in Fig. |^. 

If the two flavor QCD belongs to the same universality class as three 
dimensional 0(4) spin models, the chiral condensate should satisfy the scaling 
relation (Q) with the identification M = (^ r ^I') su b, h = 2m q a, and t = (3 — (3 ct . 
We further expect that the scaling function f(x) itself is a universal function 
15] because f(x) is determined by the universal singular structure of the 
free energy around the UV fixed point. We make a fit to the scaling function 
recently obtained for an 0(4) model |I5|1 , by adjusting (3 ct and the scales for 



t and h, with the exponents fixed to the 0(4) values. Fig. ||a shows our 
result with x^/df = 0.61. The scaling ansatz works remarkably well with 
the 0(4) exponents. The resulting (3 ct = 1.35(1) is slightly smaller than the 
value ~ 1.4 obtained by linear extrapolations of the K t line in the (/3,K) 
space (cf. Fig. []]), the on K c (cf. Fig. §), an d also the on K t |Uj. 
However, the 0(4) universality predicts |[| that (3 c {m q ) — (3 ct oc m^ 13 ^, and 
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~ (f3 c — j3 c t)' y on the K c and K t lines with 7 ~ 1.4, which implies that 
these lines should bend slightly near the chiral transition point to give a 
smaller (3 ct . Therefore, we conclude that (3 ct = 1.35(1) is consistent with the 
data. 

Finally, we test the meanfield (MF) exponents [1/(05) = 2/3, 1/5 = 1/3] 



suggested in Ref. ||16||. We find that it is more difficult to arrange a wide range 
of data on a universal curve by adjusting (3 ct . A fit with the minimum x 2 / df = 
2.5 with the MF scaling function is obtained with /3 ct = 1.54(1). However, 
the data at (3 — 1.6 are completely off the fit. Furthermore, this value of 
(3 c t is in conflict with the fact that the data of on the K c line indicates 
Pet < 1-5 (cf. Fig. §). Restricting (3 ct < 1.5, we obtain the fit shown in Fig. |b 
with x 2 /df = 3.3. In contrast with the case of the 0(4) exponents, the MF 
scaling function can not well reproduce the t/h 1 ^ 135 ^ dependence of the data 
at fixed /3, especially at small /3's. As a result, the whole data are more 
scattered than the case of the 0(4) exponents. When we decrease (3 ct less 
than 1.4, which seems natural from other data (cf. Figs. [I], §), ( 1 l , \I/) su b//i 1 ^ 5 
does not fall on an approximate universal curve. Therefore, we conclude that 
the data do not favor the MF scaling. 

The success of this scaling test with the 0(4) exponents strongly suggests 
that the chiral transition is of second order in the continuum limit. It also 
indicates that, with the improved gauge action, the chiral violation due to the 
Wilson term is sufficiently weaker than that introduced by the non-vanishing 
m q at least for m g 's at /?'s studied here. To strengthen the conclusion, a 
direct extraction of each critical exponent would be desirable. 

We are grateful to Doug Toussaint for very useful discussions and sending 
us the data of the scaling function for an 0(4) spin model. The simulations 
are performed with Fujitsu VPP500/30 at the University of Tsukuba. This 
work is in part supported by the Grant-in-Aid of Ministry of Education, 
Science and Culture (Nos.07NP0401, 07640375 and 07640376). 
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Figure V. Phase diagram for Np = 2 QCD with the RG improved gauge 
action Eq.((l|) and the Wilson quark action. Lines are to guide the eyes. 
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Figure 3: on K c . 
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Figure 4: The same as Fig. for twice the quark mass. 
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Figure 6: Best fits to the scaling function with (a) 0(4) and (b) MF ex- 
ponents (see text). The plot contains all the data shown in Fig. [| within 
the range < 2m q a < 0.8 and j3 < 2.0. Solid curves are scaling functions 
obtained in an 0(4) spin model [13 and by a MF calculation, respectively. 
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